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Abstract 

We present calculations of structure functions using a renormalization scheme consistent ex- 
pansion which is leading order in both ln(l/x) and a s (Q 2 ). There is no factorization scheme 
dependence, and the "physical anomalous dimensions" of Catani naturally appear. A relationship 
between the small x forms of the inputs F 2 (x,Qq) and Fl(x,Qq) is predicted. Analysis of a very 
wide range of data for F 2 (x, Q 2 ) is performed, and a very good global fit obtained. The prediction 
for F L (x, Q 2 ) produced using this method is smaller than the usual NLO in a s (Q 2 ) predictions for 
Fl(x,Q 2 ), and different in shape. 
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1. Introduction 

The recent measurements of F2(x,Q 2 ) at HERA have provided data on a structure function 
at far lower values of x than any previous experiments, and show that there is a marked rise in 
F 2 (x,Q 2 ) at very small x down to rather low values of Q 2 [1][2]. Indeed, the most recent measure- 
ments demonstrate that the rise persists for values of Q 2 as low as 1.5 GeV 2 . These measurements 
have led to a great deal of interest in how one should best calculate structure functions. 

The particular interest in the small x region comes about because it has long been known that 
there is potential small x enhancement of the structure functions at high orders in the strong cou- 
pling constant [3] which is known as BFKL physics: i.e. the effective splitting function governing the 
growth of the gluon Green's function at small x is of the form P(x) = J2m= 

\n m -\l/x) 

[4], where the a m are such that an asymptotic growth j- 1 " * 41112 was predicted (a s = 3a s /7r). 
This implies that one needs more than the normal fixed order in a s expansion in order to describe 
physics at small x. Qualitative studies incorporating these ideas to obtain the structure functions 
(rather than just the gluon Green's function) were in reasonable qualitative agreement with early 
HERA data [5]. 

However, it was shown to be possible to obtain very good fits to the same data whilst ignoring 
these leading ln(l/x) terms. Using the Altarelli-Parisi evolution equations at next to leading 
(NLO) order, or even leading order (LO) in a s and starting with flat [6], or even valence-like 
inputs [7], predicted a steep (though not powerlike) rise in the small x structure functions at 
scales reasonably far above the starting scale Qq, and gave a good fit as long as Qq was chosen 
to be low. This countered the BFKL approach which, after all, was derived using a less well- 
defined theoretical framework than the renormalization group approach based on the factorization of 
collinear singularities and ignored everything but the leading ln(l/x) terms. An approach somewhat 
intermediate between these extremes is also used, i.e. the fixed order evolution beginning from 
inputs for the parton distributions of the form x _1_A at small x, with A ~ 0.2, partially justifying 
the relatively steep input (significantly steeper than that expected from non-perturbative physics, 
if somewhat less than a s 41n2) from BFKL physics, e.g. [8]. 

Recently, due to the work of Catani and Hautmann [9], it is possible to use the leading in 
ln(l/x) expressions for the anomalous dimensions and coefficient functions within the well-defined 
renormalization group approach. Using the /c^-factorization theorem [10] [11] [9] they verified the 
form of the gluon anomalous dimensions 7° 9 (iV, a s ) and 7* (iV, a s ) (where j(N) = J x N P{x)dx) 
and derived expressions for jj^N, a s (Q 2 )) and r yj g (N,a s (Q 2 )) in certain factorization schemes 
(since 7//( s )(^ a s(Q 2 )) = the quark anomalous dimensions are a power of iV -1 down on those 
of the gluon). They also derived expressions for the coefficient functions C£ 1 (A r , Q 2 ), C[ 1 (A r , Q 2 ), 
Cf 1 (AT, Q 2 ) and C{ 1 (A r , Q 2 ) (all zeroth order quantities being zero except C^q, which is unity). 
This facilitated calculations of structure functions within the normal renormalization group frame- 
work whilst including the leading ln(l/x) terms. Calculations were performed using these terms 
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[12]-[14] and comparisons with data made. The calculations used different methods of solution, 
made rather different assumptions and used (different) ansatze for unknown terms. Consequently, 
differing results were obtained. The conclusions which could be drawn regarding the inclusion of 
the leading ln(l/x) terms depended very much on which approach was taken. However, including 
these terms did not seem to improve the best fits for the small x data using one- or two-loop 
evolution [12] [14]. Indeed, many ways of inclusion significantly worsened the fits, particularly if 
they were more global, i.e. constrained by large x data [12]. Also, there seemed to be a very strong 
dependence on the factorization scheme used in the calculations when including the leading ln(l/x) 
terms [12] [14] [15]. 

The high precision of the most recent HERA data constrains theory far more than previously. 
The best recent global fits seem to come from those intermediate approaches which use NLO 
perturbation theory with a quite steep (unexplained) input for the singlet quark with A ~ 0.2 and 
a similar form of small x input for the gluon [16] (unless Qq is less than ~ 4GeV 2 , in which case 
it must be flatter or even valence-like). Fixed order perturbation theory with flat or valence-like 
inputs and low Q$ fails at the lowest x values, and for fits to the small x data alone relatively 
steep inputs for the singlet quark, i.e. A <j 0.2, seem to be required [17]. Approaches including the 
resummed terms now seem to fail [18] in practically all factorization schemes. 

In this paper we will demonstrate that the apparent failure of approaches using the leading 
ln(l/x) terms, and certainly the factorization scheme dependence, is due to incorrect methods of 
incorporating these terms. The correct leading order renormalization scheme consistent (RSC) 
calculation naturally includes leading ln(l/x) terms in a form which has already been derived by 
Catani and called "physical anomalous dimensions" [19]. It also provides some limited predictive 
power at small x. We will discuss this method of calculation, then make detailed comparisons to 
data, and with the aid of the new HERA data demonstrate that this calculation leads to a very 
good global fit to all F 2 (x,Q 2 ) data. Indeed, the complete RSC calculation, including leading 
ln(l/x) terms, is clearly preferred by the latest data. 

2. The Renormalization Scheme Consistent Expansion. 

For simplicity we work in moment space in this paper, i.e. define the moment space structure 
functions by the Mellin transformation, 



Jo 

and similarly for the coefficient function. As with the anomalous dimension we define the moment 
space parton distribution as the Mellin transformation of a rescaled parton density i.e 




(2.1) 




(2.2) 
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The moment space expression for a structure function is then 



F(N, Q 2 ) = VJ C a (N, a s (Q 2 ))f a (N, Q 2 ), (2.3) 

a 

and the parton distributions evolve according to the perturbative renormalization group equation 

^^ = £7afcMQ 2 ))iUQ 2 ), (2.4) 

b 

where we choose both the factorization and renormalization scale to be equal to Q 2 . The coefficient 
functions and anomalous dimensions are factorization scheme dependent, of course. 

There are two independent structure functions F 2 (N,Q 2 ) and Fl(N,Q 2 ) which have both 
singlet and nonsinglet contributions. In general we may write 

/ N f 

Fi(N, Q 2 ) = — ( J] e 2 )F?(N, Q 2 ) + F t NS (N, Q 2 ) (i = 2, L), (2.5) 

/ \ j=1 J 

where in terms of coefficient functions and parton densities 

if (AT, Q 2 ) = C{ (N, a s )f s (N, Q 2 ) + Cf(N, a s )g(N, Q 2 ), 

N f 

if s (iV, Q 2 ) = Cr(N, a.) VJ e 2 f q N 3 S (N, Q 2 ). 



(2.6) 



Nf is the number of active quark flavours (we only consider massless quarks), and f (N,Q 2 ) and 
fq. s (N,Q 2 ) are the singlet and nonsinglet quark distribution functions respectively. 

In order to devise an expansion scheme for the calculation of these structure functions which 
is useful at both large and small x, we would a priori expect that we would need to use the known 
anomalous dimensions and coefficient functions at low orders in both a s and in the leading ln(l/x) 
expansion. There have already been some methods along these lines; however, these have suffered 
from scheme dependence. This is clearly incorrect since we do not expect factorization scheme 
dependence in a well ordered calculation of a physical quantity. 

As already mentioned, Catani has shown how one may obtain factorization scheme independent 
results, even at small x, by the use of factorization scheme independent, or physical, anomalous 
dimensions. We refer to his papers [19], or for a slightly different presentation [20], for details. 
Very briefly, using (2.3) one writes parton distributions in terms of structure functions and coeffi- 
cient functions, and substituting into (2.4) leads to evolution equations for the structure functions 
themselves in terms of physical anomalous dimensions T 2 2(N,a s ), T 2 L{N,a s ), etc. These can be 
expressed in terms of coefficient functions and anomalous dimensions in any factorization scheme, 
but in scheme independent combinations. In order to perform our calculation of structure functions 
we do not work with these physical anomalous dimensions from the outset. Rather we will simply 
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demand a complete leading order, renormalization scheme consistent (RSC) calculation. This leads 
to expressions which are unique, up to nonperturbative inputs, and which naturally contain the 
physical anomalous dimensions. 

To begin, let us consider what we normally mean by "consistency with renormalization scheme 
dependence". In the loop expansion the order of a term is determined simply by its order in 
ct s , and in the leading ln(l/x) expansion (ln(l/x)) _1 is put on an equal footing to a s . In both 
forms of expansion one demands that, once we choose a particular renormalization scheme and 
work to a particular order in this renormalization scheme, we include all terms in our expressions 
for the structure function which are of lower order than the uncertainty due to the freedom of 
choice of renormalization scheme (i.e. the uncertainty in the definition of the coupling constant), 
and no others. 1 If working with the n-loop coupling constant, the ambiguity in its defintion due 
to renormalization scheme uncertainty is of order a" +1 . Thus, the uncertainty when working 
to n-th order is the change in the leading order expression under the change of coupling a s — » 
a s (l + ecu"). Hence, the uncertainty in the whole structure function is of the order of the change 
of the leading order part under such a change in the coupling, and the n-th order renormalization 
scheme independent expression includes all complete terms of lower order than this change. 

This definition gives a well-defined way of building up a solution to the structure functions, 
but relies upon the definition of a given expansion scheme. It leaves an ambiguity about how we 
define leading order expressions and about how we define the order of terms compared to this 
leading order. The loop expansion and leading ln(l/x) expansion are just the two most commonly 
used examples. Both have potential problems: in the former one does not worry about the fact 
that the large ln(l/x) terms can cause enhancement of terms which are higher order in a s at small 
x, and in the latter one does not worry about the fact that at large x, especially as Q 2 increases, 
it is the terms that are of lowest order in a s which are most important. Hence, one would think 
that both have limited regions of validity. 

The shortcomings of these two expansion methods come about because even though a given 
order contains no terms which are inconsistent with working to this order in a particular renormal- 
ization scheme, in neither does it include every one of the terms which are consistent with working 
to a given order in the renormalization scheme. In each expansion scheme some of the terms ap- 
pearing at what we call higher orders are not actually subleading in a s to any terms which have 
already appeared. Thus, although (for a given expansion method) these terms are formally of the 
same order as uncertainties due to the choice of renormalization scheme, they are not terms which 
are actually generated by changes in renormalization scheme. 2 

1 Following this prescription one automatically obtains factorization scheme independent expressions in 
both the above expansion schemes. This is well known in the loop expansion, see [20] for a discussion of 
this and the more complicated case of the leading ln(l/x) expansion. 

2 Similarly, they are not generated by changes in renormalization scale. 
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In order to demonstrate this point more clearly we consider a simple toy model. Let us imagine 
some hypothetical physical quantity which can be expressed in the form 

oo oo oo oo , ( 1 

H(N,a s (Q 2 )) = <*s(Q 2 ) £ a mn N n J>* S (Q 2 ) £ ) " (2 " 7) 

m=l n= — m i=0 j = l — i ^ ' 

The first way of writing H(N, a s (Q 2 )) as a power series corresponds to the loop expansion, where 
we work order by order in m, out to m = k, and use the /c-loop coupling. The second corresponds to 
the leading ln(l/x) expansion, where we work order by order in i, out to i = I, and use the (I + 1)- 
loop coupling. Let us, for a moment, consider the leading order expression in the loop expansion, 
a s{Q 2 ) J2^=-i a inN n . The coupling is uncertain by 0(a 2 (Q 2 )) and hence the uncertainty of the 
leading order expression (i.e. the change due to a change of the coupling) is ~ a 2 {Q 2 ) Y^=-i bi n N n . 
We see that there is no change of any sort with a power of N less than —1, and hence any such term 
is not really subleading. Similarly, the uncertainty of the leading order expression in the leading 
ln(l/x) expansion contains no terms at first order in a s (or with positive powers of 7V), and such 
terms are not really subleading either. The full set of terms in the combination of both leading 
order expressions is genuinely leading order, and renormalization scheme independent by definition. 

Perhaps the best way in which to write our expression for H(N, a s (Q 2 )) in order to appreciate 
these points is 

oo oo oo oo 

H(N,a s (Q 2 )) = £ iV m J> mn <(Q 2 ) + £ N~ m £ c mn <(Q 2 ), (2.8) 

m= — 1 n=l m=2 n=m 

i.e. as an infinite number of power series in a s (Q 2 ), one for each power of N. Each of these series in 
a s (Q 2 ) is independent of the others, and the lowest order in a s (Q 2 ) of each is therefore renormal- 
ization scheme independent and part of the complete leading order expression for H(N,a s (Q 2 )). 
The full leading order expression for H(N,a s (Q 2 )) is therefore 

oo oo 

H (N,a s (Q 2 ))= £ N m c ml a s (Q 2 ) + Y J Cm m N- m aT(Q 2 ) 

(2.9) 



OO OO r 

^a s {Q 2 ) aonN n + Y,b 0j ( 

n=-l 1=2 ^ 



a s (Q 
N 



Hence, the combined set of terms considered to be leading order in both the previous expansion 
schemes comprise the full set of renormalization scheme invariant, and thus truely leading order, 
terms. By considering H(N,a s (Q 2 )) written as (2.8), and considering a redefinition of the coupling 
constant, a s (Q 2 ) — > a s (Q 2 ) + (D(a™(Q 2 )), we see that the n-th order expression for H(N,a s (Q 2 )), 
which should be used with the n-loop coupling constant, is the sum of the first n terms in each of 
the power series in a s (Q 2 ). 
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Similar arguments have already been applied to the anomalous dimensions and coefficient 
functions ([10] and particularly [12]). Here we take a strong viewpoint and insist that the complete 
renormalization scheme consistent expressions, with no artificial supression of leading ln(l/x) terms, 
must be used. Futhermore, the expressions must be those for the physical structure functions, not 
for the factorization and renormalization scheme dependent coefficient functions and anomalous 
dimensions. 

When considering the real structure functions the situation is a great deal more complicated 
technically than the toy model, but the principle is the same. One complication is that the structure 
functions are combinations of perturbative evolution parts and input parts (which are perturbative 
with nonperturbative factors), rather than one simple power series in a s (Q 2 ). However, the physical 
consequence of the factorization theorem is that there is some fixed nonperturbative factor for 
the structure functions which we cannot calculate in perturbation theory, but we can predict the 
perturbative evolution of the structure functions in terms of this factor. Thus, we choose two 
independent variables for each structure function, the input at some starting scale Qq, and the 
evolution away from this starting scale. We calculate the lowest order RSC expression for each of 
these, and combine to give the full LO expression. In this paper we only do this for the singlet 
structure function, since this is hugely dominant at small x. The procedure is also complicated 
here because the evolution of Fj? and is coupled, but the expression for each is a sum of terms 
consisting of input and evolution parts, and for each term we take the LO expression for the input 
and for the evolution. 

Using the above prescription it is relatively straightforward, but rather involved to calculate 
the full leading order RSC expressions. The full details appear in [20], here we just present the 
results. In order to do this most succinctly we express the result in terms of physical quantities. 
Hence, in order to explain notation (which is similar to that in [19]), and also slightly elucidate the 
form of the final expressions we consider the solution to the 1-loop renormalization group equation 
first (this being part of the full solutions). Using boundary conditions F®' 1 (N,Qq) = Fl(N) and 
F° ' l (N, Qq) = F 2 (N), where Fi(N) are nonperturbative inputs, and the superscript 0,/ denotes 
one-loop quantities, we may write the solution for the longitudinal structure function as 

Fr(N,Q*) = F^(N)(^) +P L l '^ N \^)) • ( 2 - 10 ) 

where F L (N, Q 2 ) = F L (N, Q 2 ) / (a s {Q 2 ) /2ir), f°.'.+(-)(JV) are the two eigenvalues of the 0(a s ) 
physical anomalous dimension matrix, which are the same as those of the 0(a s ) parton anomalous 
dimension matrix, divided by boa s (Q 2 ). (The superscript S is dropped for the rest of this section.) 
Having chosen to write the solution for F®' l (N,Q 2 ) in this way we may then write F°' l (N,Q 2 ) as 

F^\N,Q 2 ) = e^N)F^(N)[^ j j + e~ (N)F° L > 1 >- (N) (2.11) 
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where e+ (N)F®' 1 ' + (N) + e~(N)F°' 1 ' (TV) = F 2 (N). The e+^(N) come from the eigenvectors 
of the 0(a s ) parton anomalous dimension and the 0(a s ) longitudinal coefficient functions. In 
practice 



F o,i,+ {N) = f l[n) _ 36 J Nf F 2 (N) + 0{N), F°/>-{N) = - n ^' J F 2 (N) + 0{N), 



36 - 8Nf 



27 



27 



F o,i,+ {N) 



N 



F L (N) 



36 - 8N f 
27 

A/.+C-), 



F 2 (N) + 0(iV 2 ), F%' l '~(N) = F 2 (N) + 0(iV), 



(2.12) 



where F 2 °' i ' +( " ) (N) = (N)F^ ( ~' {N). 

Taking into account the leading ln(l/x) terms as well, the expressions acquire additional terms. 
Explicitly we obtain 



L,RSC,0 



(N,Q 



2, ._ <x 8 (Ql) 



2vr 



^M" f °' i,+(JV) - 1 



exp 



36 - 8N 



27 



for the longitudinal structure function, and 
dF 2 (N,Q 2 



« S (Q 2 ) b a 2 (q 2 
F 2 (N))(exp[ln(Q 2 /A LL )r LL (N,a s (Q 2 ))] - 1) 



(2.13) 



cUnQ 2 



= a s (Q 2 ) 



RSC,0 



e~ (N)T°' l '~ (N)F^' l '~ (N) 



0,1, 



a s (Q 2 ) 



+ (e^N)T^(N)F^(N) -T^ L (N)(p L (N) - ^-J^F, 



(N) 



+ Tl L (N,a s (Q 2 ))(p L (N) - ( 36 2 * Nf )F 2 (N)\ exp[ln(Q 2 M LL )r° L (iV, a s (Q 2 ))] 



x 



exp 

and 

F 2 ,rsc,o(N,Q 2 ) =F 2 (N) 



«.(0g) r| L (JV,Q 8 (g 2 )) 



6 a 2 (<? 2 ) 



da s (g 2 ) 



a s (Q 2 ) 



(2.14) 



+ a - (9 ° ) iSSii ( #L(iV) " (36 2 j Nf) ^ N ))^MHQl/A LL )rUN,aM))] - 1) 



+ li l (^r,/.lLL)n.(^)(c-(.V)r (l -' ^(N^+iN.Qi) ->re-(N)T^-(N)F^-(N,Q 2 ) 

'36 - 8iV/ 



rl>° L (N)[F L (N) 



27 



F 2 (iV) 



(2.15) 



for F 2 (N,Q 2 ). 



We should explain the terms in these expressions. T LL (N, a s ) is the gluon anomalous dimension 
at leading order in ln(l/x), which it turns out governs the small x evolution of Fl(N, Q 2 ), as seen 
in (2.13). It is also identical to the physical anomalous dimension, hence the notation. T LL (N,a s ) 
is T LL (N, a a ) minus its one-loop component, which appears in T°' l ' + (N). a s Tl L (N, a s ) is the 
leading order in ln(l/x) term governing the evolution of F 2 (N,Q 2 ) in terms of Fl(N,Q 2 ). It is 
a power of a s up on T° LL (N,a s ), but is sufficiently leading (i.e. is not subleading to any other 
contributions) to make an appearance in the input terms in (2.14) and (2.15). It is given by 

jlJN, a s ) + 7° C\ JN, a s ) 

and thus is equal again to one of the physical anomalous dimensions in [19]. (T^l (N) is the one- 
loop contribution to Tl, L (N, a s ), and must be subtracted in some places in (2.14) and (2.15) in 
order to avoid double counting.) However, once again we stress that these anomalous dimensions 
are not needed to derive these expressions, but that they naturally appear in the end results. 

The term exp[ln(Ql/ALL)T° L L (N, a s (Q 2 ) ))] appears in some of the input parts in our expres- 
sions. This leading ln(l/x) contribution to the inputs is derived by demanding that the form of 
our expressions is invariant under changes in the arbitrary starting scale Qq at the order at which 
they are calculated, i.e. the expression for the structure functions as a whole is genuinely of leading 
order. It is easy to see that the variation of this input term under a change in Qq cancels the leading 
order change in the leading ln(l/x) evolution. This procedure leaves us with an unknown scale 
All, but at this scale the inputs become the nonpertubative inputs, and we would therefore expect 
All to be the sort of scale where perturbation theory starts to break down. The other terms in 
(2.15) proportional to Io.(Qq/All) are likewise required to make the full expression invariant under 
changes in Qq at leading order. 

Having obtained the full leading order RSC expressions for (dF 2 (N, Q 2 )/d In Q 2 ) and 
F 2 (N,Q 2 I ) we integrate (dF 2 (N,Q 2 )/dlnQ 2 ) RS c,o from Qq to Q 2 and add to the input 
F 2 ,rsc,o(N,Qq) in order to get our expression for F 2 (N,Q 2 ). This is essentially because it is 
the derivative of F 2 (N,Q 2 ) that begins at first order in a s and thus is a truly perturbative quan- 
tity. The difference between using this prescription and using the leading order F 2 (N,Q 2 ) directly 
is small. 

It is easy to check that under a change in the coupling, a s — > a s + 5a 2 the change of each of 
our expressions is of higher order in a s than any terms appearing. Thus, we have full leading order, 
including leading ln(l/x) terms, renormalization scheme consistent expressions for the structure 
functions. These are significantly different from both the one-loop expressions and the leading 
ln(l/x) expressions, although they clearly reduce to them in the appropriate limits. All terms in the 
expressions (2.13)-(2.15) are renormalization scheme as well as factorization scheme independent, 



as we would hope, and they are the appropriate full expressions to use with the one-loop coupling 
constant. 

Finally, we notice that this method of solution leads to a certain amount of predictive power. 
We know the precise form of the structure function inputs in terms of the nonperturbative inputs 
(which we imagine should be quite flat at small x). Thus, up to the absolute normalization and 
the scale All, we have predictions for the small x form of the inputs for Fl(N, Q 2 ) and F 2 (N, Q 2 ) 
(as well as for (dF 2 (N, Q 2 ) / din Q 2 )) . The normalization is fairly well set by the large x data, 
and we would expect All ^ lGeV 2 , so there is an estimate for the small x form of each input. 
Moreover, the unknown elements are the same for each input, so there is a strong prediction for 
the relationships between the small x inputs. However, the scale Qq which should be chosen is 
not determined. Nevertheless, it is a considerable consistency requirement that the relationships 
should be true for any choice of Qq, and hopefully they can be well satisfied over a wide range of 
Qq. In order for this to be true, a s (Ql) cannot be too sensitive to Qq, so we would not expect Qq 
to be particularly low. Also, since we can largely choose our structure functions at Qq, but then 
have no freedom in how we evolve up and down in Q 2 , we would imagine that when performing a 
fit it is advantageous if Qq is near the centre of the range of our data. 



3. Fits to The Data. 

We use the expressions (2.13)-(2.15) to calculate the x-space singlet structure functions. The 
nonsinglet structure functions are calculated using the normal one-loop prescription. By combining 
the singlet and nonsinglet components and varying the free parameters (Qq, All, and the soft 
inputs for Fl(x,Q 2 ) and F2{x,Q 2 )), we obtain the best fit for the available F 2 structure function 
data. 3 We note that the input F 2 (x, Qq) and the evolution dF 2 (x, Q 2 )/d\n Q 2 are forced by (2.14) 
and (2.15) to be trivially related at small x, which is not the case when working at fixed order in 
a s . The one-loop value for A Nf=4 is chosen to be lOOMeV. This precise value is not determined 
by a best fit, but a value near this is certainly favoured. The published values of F 2 (x,Q 2 ) are 
altered to take account of the fact that our predictions for Fl(x, Q 2 ) are not the same as (i.e. are 
somewhat lower than) those used by HI and ZEUS in their determination of F 2 (x, Q 2 ). Thus, the 
F 2 (x, Q 2 ) values are a little lower for the largest values of Q 2 jx than in [1] and [2]. 

3 In practice a variant of the program used by MRS is used for the fit, and inputs for the gluon and 
quarks of the standard form are specified. To calculate structure functions we use an effective factorization 
scheme which is of DIS type for F2(x,Q 2 ), where the longitudinal coefficient functions are the one-loop 
expressions, and the resummed anomalous dimensions and small x inputs are chosen to produce expressions 
for the structure functions matching (2.13)-(2.15) as closely as possible. For simple inputs where exact 
analytic expressions for the x-space forms of (2.13)-(2.15) can be found checks are made with the results 
of the evolution program, and any discrepancies are always much smaller than the errors in the data. 
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We treat the heavy quark thresholds rather naively. Taking m 2 = 3GeV 2 and to 2 . = 20GeV 2 , 
we simply change the number of active quark flavours discontinuously at these values of Q 2 : since 
F[(N) has a large component proportional to Nf, Fl(x,Q 2 ) is discontinuous at these values of 
Q 2 , as is dF 2 (x, Q 2 )/d\n(Q 2 ). a s (Q 2 ) is continuous at the thresholds, being defined by 

a s , n {Q 2 ) = « s ,n+i(Q 2 ) (l + as ^ Q2) ln(m 2 +1 /Q 2 )) . (3.1) 

This treatment of quark thresholds is consistent with the decoupling theorem, in so much that it 
guarantees the correct expressions far above or below thresholds [21]. It is clearly unsatisfactory 
near the thresholds and must be improved. However, the prescription has little effect on F 2 (x, Q 2 ) 
in the region of the fit: in the curves for F 2 (x, Q 2 ) the 6-quark threshold is barely noticeable, while 
the kink at the c-quark threshold only really affects a handful of data points at very small x, tending 
to hinder the fit (see fig. 1). 

The result of the best fit using the leading order, including leading ln(l/x) terms, RSC expres- 
sions (henceforth refered to as LO(x)) with Qq = 40GeV 2 is compared with fits obtained using the 
standard two-loop method, where Ri allows A^|T to be free (giving AlgT 4 = 241MeV) and R 2 
fixes A^lT 4 = 344MeV to force a better fit to the HERA data. The new NMC data [22] for F% p 
and F^ is used with a lower Q 2 cut of 2GeV 2 . 4 The results are shown in table 1 (full references 
for the experimental data can be found in [16]). As one can see, the LO(x) scheme independent 
fit is much better for the HERA data (even when compared to R2), much better for the BCDMS 
data (even when compared to Ri) and similar in standard for the rest of the data. The overall fit 
is ~ 200 better for the whole data set. The results of the fit to the small x data is shown in fig. 1. 

The leading order renormalization scheme consistent expressions clearly provide a very good 
fit to the F 2 (x,Q 2 ) data. The fit shown is for the particular starting scale Qq = 40GeV 2 , but the 
quality of the fit is extremely insensitive to changes in this scale (where we allow All to be a free 
parameter for each Qq), as we expect from the method of construction of the solutions. The fit is 
essentially unchanged over the range 20-80GeV 2 , and we choose 40GeV 2 as the (logarithmically) 
central value. When Q 2 , drops below 20GeV 2 the fit immediately gets markedly worse because of 
the discontinuity in dF 2 (x, Q 2 )/d In Q 2 at the threshold: dF 2 (x, Q 2 )/dln Q 2 suddenly becomes 
too large at small x if F^x^Qfy is the correct size to fit the data. The quality of the fit gets 
continuously worse as Qq lowers further, becoming completely uncompetitive long before reaching 
m 2 . We expect that a correct treatment of quark thresholds would lead to a smooth falling off of 
the quality of the fit, but that it would begin to deteriorate somewhere in the region of 20GeV 2 , 
due to a s {Q 2 ) ) becoming too large below this value. 

4 The MRS fits are not performed again: the % 2 for the new data is calculated using the same input 
parameters in [16] (there is little indication that these would be changed much by the new data). 
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The parameter All, which should be a scale typical of soft physics, turns out to be 0.4GeV 2 
for the fit starting at Qq = 40GeV 2 . This decreases a little as Qq increases and vice versa. For 
Qq = 40GeV 2 the soft inputs for the fit are roughly 

F[ (x) « 2.65(1 -x) 5 , Fi(x) » (1 -xf(l - 0.6x a5 + 7x), (3.2) 

where they have been forced to be flat as x —> 0. Allowing instead an asymptotic behaviour x x , 
where A & 0.08, leads to an equally good fit. The importance of the leading ln(l/x) terms can be 
judged by how they affect the fit. If, after obtaining the best fit, all terms other than those in the 
one-loop expressions are set to zero, the quality of the fit is unchanged above x = 0.3, begins to 
alter slightly below this, and is clearly much worse by the time we reach x = 0.1. Thus, the leading 
ln(l/x) terms are important by this value of x. 

It is not yet possible to extend the RSC calculation beyond the leading order due to lack of 
knowledge of NLO in ln(l/x) terms. There is hope that these will shortly become available [23], 
and when they do the NLO versions of (2.13)-(2.15) can be derived and put to use. Only then 
should the NLO coupling constant be used in any fit. As shown at leading order the ln(l/x) terms 
not present in the one-loop expressions become important above x = 0.1. However, much of this 
effect is due to the terms at C(a 2 ), so the NLO expression at fixed order in a s should be a good 
approximation to the full NLO RSC expression for x somewhat lower than 0.1, perhaps as low as 
~ 0.05 — 0.01. However, until the full renormalization scheme consistent NLO expressions become 
available, we believe that it is premature to use fits to small x structure function data to determine 
the NLO coupling constant (unless, of course, direct measurements of Fl(x,Q 2 ) and other less 
inclusive quantities at very small x turn out to verify standard two- loop predictions). The fixed 
order in a s expressions should be accurate for CCFR, BCDMS and NMC data, which after all 
are still much more precise than HERA data, and fits to these data alone will provide the best 
determination of the NLO a s (Q 2 ). 

4. Conclusion 

In this paper we have demonstrated that it is possible to derive expressions for the structure 
functions which incorporate the leading ln(l/x) terms in a way which is renormalization scheme 
consistent, and as a consequence avoids any factorization scheme dependence. We have also shown 
that these full leading order RSC expressions lead to very good fits to the data and that, futhermore, 
they are able to do so using as inputs only soft distributions for the singlet quark and gluon, i.e. all 
powerlike behaviour is generated perturbatively, and determined in terms of the nonperturbative 
flat inputs and a soft scale All- Hence, this approach provides an explanation for the form of 
the small x structure function rather than just a way of fitting it. Futher details of both the 
theory and fits, as well as other related issues, are presented in [20]. It is certainly true that the 
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calculations must be improved to take account of massive quark thresholds in a better manner, 
and work towards this end is in progress. Nevertheless, with the present treatment we feel that 
the quality of the fit and the degree of explanatory (if not predictive) power, not to mention the 
scheme independence, give strong justification for using this approach. 

However, the quality of the fit alone is certainly not such a substantial improvement on more 
standard approaches that it necessarily convinces one that this approach must be correct. In order 
to obtain some degree of verification we must obtain more experimental data. So far we have only 
probed Fl(x,Q 2 ) indirectly, i.e. it is simply related to the derivative of F 2 (x,Q 2 ) (as well as to 
the input F 2 (x,Qq) using this method). Having tied down the nonperturbative inputs and All 
and Qq from our fit to F 2 (x,Q 2 ), we have a prediction for F L (x,Q 2 ). The result of this prediction 
for the fit with Ql = 40GeV 2 is shown in fig. 2, where it is compared to the prediction using the 
NLO in a s approach and the MRSRi fit. As one can see, it is smaller than the MRSRi Fl(x, Q 2 ), 
but becomes steeper at very small x. The prediction for Fl(x,Q 2 ) is weakly dependent on the 
value of Q 2 chosen: the value at Q 2 = 5GeV 2 and x = 10 -4 varies by ±10% within our range 
of Qq (increasing with Qq), and by less than this for higher x and Q 2 . 5 The very recent results 
on F L (x,Q 2 ) for 0.01 k, x Z 0.1 from NMC [22] are matched far better by the LO(x) F L than 
the MRSRi Fl (the latter being rather large). However, it is fair to say that any problems with 
the MRSRi Fl can very probably be assigned to the treatment of the charm quark threshold, i.e. 
the predicted Fl(x,Q 2 ) in the last of [7] matches the data well. Measurements of Fl(x,Q 2 ) at 
x < 10 -2 would be a better discriminant between fixed order in a s calculations and those involving 
leading ln(l/x) terms. However, the sort of "determination" of Fl(x, Q 2 ) already performed by HI 
[25] is really only a consistency check for a particular fit, and is by no means a true measurement of 
Fl(x, Q 2 ). Real, direct measurements of Fl(x, Q 2 ) at HERA would be an important (and probably 
essential) way of determining the validity of the approach in this paper, and the genuine importance 
of leading ln(l/x) terms in structure functions. 
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The author has submitted two conference proceedings on the topic of the current paper [24] and should 
point out that these are both incomplete. In the former the expression (2.15) was not used for the input for 
F2(x,Q 2 ). Thus, the fit imposed far less constraint on Qq than the full procedure, and a value of 5GeV 2 
was used, which resulted in a prediction of Fl(x,Q 2 ) that is much too small. The latter claimed that an 
input for Fl(x, Q 2 ) a little smaller than that consistent with the full set of expressions was needed for the 
best fit, even at the optimum Qq. This was due to there being no account whatever taken of the 6-quark 
threshold. The correction makes very little difference to phenomenological results. 
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Table 1 

Comparison of quality of fits using full leading order (including ln(l/x) terms) renormalization 
scheme consistent expression, LO(x), and two-loop fits MRSRi and MRSR2. For the LO(x) fit the 
HI data is normalized by a factor of 0.995, the ZEUS data by 1.01, the BCDMS data by 0.98, the 
CCFR data by 0.95, and the rest by 1.00. 



T^vnpn m Pint 


K. 1 CI I CI 




y 2 

A 






points 


LO(x) 


Ri 


R 2 


HI F* p 


193 


128 


158 


149 


ZEUS F 2 ep 


204 


256 


326 


308 


BCDMS F% p 


174 


190 


265 


320 


NMC F£ p 


129 


124 


163 


135 


NMC F£ d 


129 


109 


134 


99 


NMC F£ n /F£ p 


85 


142 


136 


132 


E665 F 2 MP 


53 


8 


8 


8 


CCFR F% N 


66 


51 


41 


56 


CCFR F$ N 


66 


49 


51 


47 
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Figure Captions 

Fig. 1. The curves correspond to the value of the proton structure function F2(x,Q 2 ) obtained 
from the leading order renormalization scheme consistent (LO(x)) calculation at 12 values 
of x appropriate for the most recent HERA data. For clarity of display we add 0.5(12 — i) 
to the value of F 2 (x,Q 2 ) each time the value of x is decreased, where i = 1 — ► 12. The 
data are assigned to the x value which is closest to the experimental x bin (for more details 
see the similar figure displaying the two-loop fits in [16]). E665 data is also shown on the 
curves with the five largest x values. The HI and ZEUS data are normalized by 0.995 and 
1.01 respectively in order to produce the best fit. 

Fig. 2. Comparison of predictions for Fl(x, Q 2 ) using the full renormalization scheme consistent 
(LO(x)) fit and the two-loop MRSRi fit. For both sets of curves Fl(x, Q 2 ) increases with 
increasing Q 2 at the lowest x values. 
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F L (x,Q 2 ) at Q 2 = 5,10, 10 2 , 10 3 GeV 2 




